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We find a method to rewrite the equations of motion of scalar fields, generalized DBI field and 
quintessence, in the autonomous form ioxarbitrary scalar potentials. With the aid of this method, 
we explore the cosmic evolution of generalized DBI field and quintessence with the potential of 
multiple vacua. Then we find that the scalars are always frozen in the false or true vacuum in 
the end. Compared to the evolution of quintessence, the generalized DBI field has more times of 
oscillations around the vacuum of the potential. The reason for this point is that, with the increasing 
of speed 0, the friction term of generalized DBI field is greatly decreased. Thus the generalized DBI 
field acquires more times of oscillations. 
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I. INTRODUCTION 


Scalar fields play an important role in both theoretical 
physics and modern cosmology for their simple but non¬ 
trivial dynamics. In theoretical physics, they are present 
in the Jordan-Brans-Dicke theory as Jordan-Brans-Dicke 
scalar in Kaluza-Klein compactification theory as the 
radion [^; in the Standard Model of particle physics as 
the Higgs boson Q; in the low-energy limit of the su¬ 
perstring theory as the dilaton Q, tachyon @ and DBI 
(Dirac-Born-Infeld) field In cosmology, scalar fields 

are employed to model the inflaton Q, the quintessence 
dSI (for a recent review of quintessence, see Refs. [Tsj l 
and references therein), the k-essence [l^, 171, the phan¬ 
tom [I^, in order to drive the inflation of the early Uni¬ 
verse or to speed up the expansion of the late Universe. 
We note that most of the scalar fields are hypothetical 
particles. But the Higgs Boson has uniquely been discov¬ 
ered by experiments. 


Quintessence is a canonical scalar field which is as¬ 
sumed to be minimally coupled to gravity. Compared 
to other scalar fields, quintessence turns out to be the 
simplest scenario which is free of ghosts and instability 
problems. The dynamics of quintessence in the presence 
of matters has been studied in great detail for many dif¬ 
ferent potentials (l2l - [T3 . fiol - l^ . However, for a general 
quintessence potential, the equations of motion are rather 
involved. To our knowledge, one have not yet find a gen¬ 
eral method to write the equations of motion in the au¬ 
tonomous form for arbitrary potential. Thus the purpose 
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of this article is to report that we have found a way. 

The DBI field describes the dynamics of D-branes 
evolving in a higher-dimensional warped spacetime. A 
novel aspect of this field is the existence of a speed limit 
on the field space, resulting from causality restrictions 
on the motion of the branes in the bulk spacetime. The 
speed limit is enforced by the non-canonical kinetic terms 
in the DBI field. This is different from the quintessence 
whose speed Vcf) could be arbitrarily large. From this 
point of view, quintessence and DBI field are the coun¬ 
terpart of Newtonian and Special Relativity mechanics, 
respectively. The cosmic evolution of DBI field have been 
studied in Refs. [2^. These researches only apply to 
some special forms of potentials. Thus, to find a general 
method applying to arbitrary DBI potential constitutes 
the second purpose of this article. 


II. GENERALIZED DBI FIELD 


We consider a D3-brane with tension T evolving in a 
5-dimensional spacetime. The dynamics of the mobile 
D3-brane is described by the DBI action. The D3-brane 
is free to move on the internal compact Galabi-Yau man¬ 
ifold. The generalized DBI action can be written as fol¬ 
lows @ 


S = j d^x^g 

-T{i;) + V{i;)] + Sm ■ 


TW ^ 




( 1 ) 


Here T{^l)) is the warped tension of the brane and Sm 
is the action for matters localized in the D3-brane. The 
potential W{tp) could arise under the condition that the 
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brane is a non-BPS one @ or there are multiple coinci¬ 
dent branes [s^. The potential is related to the 

brane’s interactions with bulk fields or other branes. 

In order to simplify our derivations, we define the vari¬ 
able (j) as follows 


1 

7^ 






Then above action can be written as 


( 2 ) 


S = j [t {^) W (<^) v^l + 

-T (</>) + ¥(</>)] +. (3) 


constrained to be smaller than the speed of light. This is 
remarkably different from the usual quintessence which 
could have arbitrary large speed in the sense of classical 
mechanics. 

Observing Eqs. ([6]) and Eq. 0, we could absorb the 
constant k^/3 = Stt) into W,V and pm, respectively, 


w - 


(9) 

V — 

CO 

7 

(10) 

Pm ~ 

3 

Pm- ^ ■ 

(11) 


Without the loss of physical significance, we could absorb 
the term into IE(0) and V(i^), respectively. Then 
we find the action is simply 


s = w{^)./TTd;;^ + vi^) + 5 ^ 4 .) 


We shall investigate the cosmic evolution of the DBI 
field in the background of spatially flat Eriedmann- 
Robertson-Walker Universe 


ds^ = —dt^ + a [dr'^ + , (5) 

where a(t) is the cosmic scale factor. We model the mat¬ 
ter sources present in the Universe as perfect fluids. The 
perfect fluids can be baryonic matter, relativistic mat¬ 
ter and dark energy. We assume there is no interaction 
between the generalized DBI scalar field and the matter 
fields, other than by gravity. Then the Einstein equa¬ 
tions and the equation of motion of the scalar field are 
given 

3H^ = + V + pra) , 

2H + 3H^ = -K^ {-Wh - U + , (6) 


and 




W 


( 7 ) 


respectively. Here H = d/a denotes the Hubble pa¬ 
rameter and dot is the derivative with respect to cos¬ 
mic time, t. Pm and ojm are the energy density and 
the equation of state for the matter sources. We have 
ujm = — 1, 0, 1/3, -|-1 for the cosmological constant, dust 
matter, relativistic matter and stiff matter, respectively. 
In this paper, we shall consider the case of dust matter, 
uJm = 0. 1/0 and denote the derivative with respect 
to (j). 7 is defined by 


7 = , , (8) 

which has the physical meaning of the generalized 
Lorentz factor. It is apparent the speed of scalar cj) is 


Then above equations of motion turns out to be 


H^ = Wj + V + pm. 

2H/3 + H^ = W/-f + V - UJmPm , 
In..; I Va I W a, 


= 0 


( 12 ) 

(13) 

(14) 


Given the scalar potential V{(j)) and the equation 

of state Wm, we are left with three variables, a(t), pm and 
the DBI field 4>(t). Then we have three variables and 
three independent differential equations. So the system 
of equations is closed. 

It is rather difficult to find the analytic solutions to 
the equations of motion (12-14). Hence in order to solve 
them numerically, we had better rewrite them in the au¬ 
tonomous form. To this end, we introduce the following 
dimensionless quantities 



U = 

wt ’ 

(15) 

H ' 

Q = ^ 
lUi ’ 

(16) 

z-^ 

c _ 1/0 
ul ’ 

(17) 

fV = Ina . 


(18) 


We see U, Q, S are the functions of DBI field, (j). So 
they can be expressed as the function of YjZ. Now we 
have only three variables, namely, X, T, Z and it is 
sufficient for us to derive the corresponding three inde¬ 
pendent differential equations. We remember that the 
system of equations of X, Y Z are an autonomous sys¬ 
tem if and only if C/, Qj S are expressed as the functions 
of X, Y, Z. For simplicity, we assume 


_U 

W 


(19) 


We emphasize that one could in general assume VjW = 
with F{(j)) other simple functions, for example, 
F{(j)) = e°“^, \n(j) and so on. Then we could obtain 

cj) = (j) iyjWj = (j) jZjYj from this equation. Similarly, 
U, Q, S could be expressed as the function of ZjY. 
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Now with the aid of this assumption, we are able to 
write the Eqs. (12-14) in the autonomous form and deal 
with any scalar potential, V{if) in the calculations. In 
this article, we shall focus on the scalar potential V{(l)) 
with the expression of 


V = Vo +-^ ((/> - ai) ((/) - 02) (c/> - as) 

■ ((/) - 04 ) ((/> - as) {(j) - ae) , ( 20 ) 


where Vi and Oi are all positive constants. The reason for 
this choice of potential is that what we want to study is 
a potential with multiple vacua. Scalar fields with mul¬ 
tiple vacua are very interesting because they have both 
theoretical origin in string landscape an d theoretical 
study in Coleman-De Luccia tunnelling [32j| . 

We note that the case of the well-known AdS throat, 
W{(j)) = has been included in the desired one. As 
an example, we put Vq = 2 , Vi = 10 ^, ai = 1 , 02 = 
11/10, as = 9/8, 04 = 39/25, 05 = 157/100, og = 2 in 
the following discussions. 

In Fig.lU we plot the potential V{(j)) with respect to />. 
There are two local maximum (.^1 = 1.11, ^2 = 1.56), 
one real vacuum (173 = 1.863) and two false vacuum 
(cti = 1.029, (72 = 1.305) for the potential. Physically, 
the scalar field would roll down the potential and then 
passes through the first [cri) and the second {(T 2 ) false 
vacuum. Finally, it arrives at the real vacuum (< 73 ). The 
detail of the trajectory is closely related to the initial 
velocity, (j)i (with the initial value, (pi fixed). When the 
initial speed pi is small enough, the scalar would acquire 
damped oscillations due to the Hubble friction in the first 
false vacuum and finally it is frozen in this vacuum. How¬ 
ever, with the increasing of initial speed, the scalar could 
cross the first local maximum (^ 1 ) and finally is frozen 
in the second false vacuum (( 72 ). With even much larger 
initial velocity, the scalar could cross the two local max¬ 
imum (^1 and ^ 2 ) and Hnally dwells on the real vacuum 
( 0 - 3 ). Since the speed of the scalar is constrained to be 
smaller than the speed of light, the scalar can not climb 
the hill with arbitrary height. Due to the Hubble fric¬ 
tion, we expect the fate of the scalar is to dwell on the 
real vacuum. In what follows, we shall show theses points 
numerically. 

Using the dimensionless variables dehned in (15-18), 
the equations of motion (12-14) can be written in the 
following autonomous form 


dX 

dN 

dY 

dN 

dZ 

dN 


-\-3X-^-UY-^-QY 

ryZ, /-yO ryZ. 

\QXY^-Y§2( 


( 21 ) 



FIG. 1: The DBI potential V{p) with respect to p. There are 
two local maximum (U = 1.11, ^2 = 1.56), one real vacuum 
( 0-3 = 1.86) and two false vacuum (ci = 1.03, (72 = 1.31) for 
the potential. Physically, the scalar field would roll down the 
potential and then passes through the first (cri) and the sec¬ 
ond ( 72 ) false vacuum. Finally, it arrives at the real vacuum 
(73). 


where U, Q, S are the functions of Y/Z and 
1 


7 = 


Vi-a:2 ’ 


( 22 ) 


H _ o 1 ~ 7 ^ 3 


= 3F^ • 

27 


- 2 (1 - 7^^" - Z^) . (23) 


The Friedmann equation becomes the constraint equa¬ 
tion 


1 = 




Pm 

IP ■ 


Vi- a:2 

The equation of state w of the DBI scalar field is 


-1/7- z2/r2 

W = 




(24) 


(25) 


In Table I, we present the properties of the five fixed 
points for the scalar field. The points (a, c, e) corre¬ 
spond to the false vacua (cti, 72 ) and real vacuum ( 73 ). 
The three points are stable spirals. On these epoches, 
the scalar field behaves as a damping oscillator with the 
equation of state of firstly behaving as the dust mat¬ 
ter and then oscillating approaching —1 (see Fig. (5)). 
The points (b, d,) correspond to the two local maximum 
(Cl, C 2 ) and they are saddle points. On these points, the 
speed of the scalar field exactly vanishes and the DBI field 
acquires the equation of state of cosmological constant. 

In Fig. (2-4), we plot the evolution of the speed, p of 
the generalized DBI field with p. We fix the initial values 
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1 Name 

|A 

y 1 ^ 1 (j) 1 

Stability [ 


\w^ \ 

1 (a) 

1 0 

0.70 1 0.72 1 (Ti = 1.03[Stable spiral (attractor) j 

1 

i-li 

1 (b) 

1 0 

0.67|0.74|G = 1-111 

Saddle point [ 

1 

i-li 

1 (c) 

1 d 

0.6110.791 (72 = 1.31 [Stable spiral (attractor) [ 

1 

i-li 

1 (d) 

1 0 

0.54 [0.84 [^2 = 1.56 [ 

Saddle point [ 

1 

i-li 

1 (e) 

0 

0.47[0.88[(73 = 1.86[Stable spiral (attractor) [ 

1 

-1 


TABLE I: Properties of the critical points for the scalar potential given by Eq. (20). 


of (j), (j) and Pm !By this way, the initial values of Y 
and Z are determined. Fig. (2) shows that when the 
initial speed (f)i is small enough, the scalar would acquire 
damped oscillating due to the Hubble friction in the false 
vacuum (cri) and finally it is frozen in this vacuum. With 
the increasing of initial speed, the scalar crosses the first 
local maximum (^i) and finally is frozen in the second 
false vacuum ( 172 ) (see Fig. (3)). With even much larger 
initial velocity, the scalar crosses the two local maximum 
(^1 and ^ 2 ) and finally dwells on the real vacuum ((T 3 ) 
(see Fig. (4)). 

In Fig. (5), we plot the evolution of the equation of 
state w of the generalized DBI field corresponding to 
Fig. (4). We see the generalized DBI field previously be¬ 
haves as the dust matter and then oscillating approaches 

— 1 after a sufficient long time. The reason for oscillating 
of w can be understood as follows. Eq. (25) tells us when 
the speed of DBI field vanishes, the equation of state is 

— 1. Fig. (4) shows there are many times for the vanish¬ 
ing of speed during the damped oscillating. Every time 
the generalized DBI field acquires vanishing velocity, its 
equation of state is — 1 . 


III. QUINTESSENCE 

In this section, we shall present the method for dealing 
with arbitrary quintessence potential. As an example, 
we would explore the cosmic evolution of quintessence 
field with multiple vacua. To this end, let’s focus on the 
potential as follows 

V = (26) 

where Vg, 6 , are positive constants. As an example, 
we consider, Vq = 1, 5=1, oi = 1/2, 02 = 3/5, 03 = 
1, 04 = 17/10, 05 = 2, 06 = 5/2. 

In Fig. [51 we plot the potential U(^) with respect to </>. 
There are two local maximum (,^1 = 0.813, ^2 = 1-862), 
one real vacuum (173 = 2.342) and two false vacuum 
(cti = 0.545, (72 = 1.355) for the potential. Physically, 
the scalar field would roll down the potential and then 
damped oscillates between these vacua. Given an initial 
finite velocity (/ and field value (/o, the fate of the scalar 
is expected to dwell on the one of the vacuum due to the 
Hubble friction. In what follows, we shall show theses 
points numerically. 



FIG. 2: Evolution of the speed (X = (/) of generalized DBI 
field with (j). The point {cj) = a\ = 1.03, A = 0) is a stable 
spiral and thus an attractor. In this case, the generalized DBI 
field behaves as a damping oscillator in the false vacuum (cti) 
and finally is frozen. 


The Einstein equations and the equation of motion of 
the quintessence are given by 

(i 02 + 1 /+ , 

2H + = —K^ ^ 2 ^^ — V Y oJmPm^ , (27) 

and 

-b + V^=0, (28) 

respectively. Here pm and uim are the energy density 
and the equation of state for the matter sources. In this 
section, we also consider the case of dust matter, ujm = 0 . 

Observing Eqs. (l27l) and Eq. (l28ll . we could absorb the 
constant k?/3 (k^ = Stt) into and pm, respectively. 



















5 



FIG. 3: Evolution of the speed {X = (j>) of generalized DBI 
field with (j). The point (0 = <72 = 1.31, X = 0) is a stable 
spiral and thus an attractor. In this case, the generalized DBI 
field behaves as a damping oscillator in the false vacuum ((T2) 
and finally is frozen. 



FIG. 5: The evolution of the equation of state for the general¬ 
ized DBI scalar held. It behaves as the dust matter at higher 
redshifts and oscillating approaches —1 at the lower redshifts 


Then above equations of motion turns out to be 



FIG. 4: Evolution of the speed {X = <j>) of generalized DBI 
held with (j). The point (0 = <72 = 1.86, X = 0) is a stable 
spiral and thus an attractor. In this case, the generalized DBI 
held behaves as a damping oscillator in the real vacuum (<73) 
and hnally is frozen. 


H^ = ^-^ + V + pm, (32) 

2H/3 + H^ = -^^ + V - WmPm , (33) 

<(i + 3if<^+il<0 = O. (34) 

Given the scalar potential I^((/>) and the equation of state 
ujm, we are left with three variables, a(t),pm and the 
quintessence (/)(t). Then we have three variables and 
three independent differential equations. So the system 
of equations is closed. 

The same as the DBI case, it is rather involved to find 
the analytic solutions to the equations of motion (32-34). 
Hence in order to solve them numerically, we had better 
rewrite them in the autonomous form. To our knowl¬ 
edge, one have only explored some special form of the 
quintessence potential, namely, the exponential potential 
[3 [13 1 power-law type potential 0113 , the Albrecht 
and Skordis potential |33j an so on. For arbitrary poten¬ 
tial, one have not yet find a general method to write the 
equations of motion in the autonomous form. In what 
follows, we shall propose a method that can be used to 
deal with arbitrary potentials. 

To this end, we introduce the following dimensionless 
quantities 


02 - 

-^^2 6 

(29) 

V — 


(30) 


k 2 

Pm ~ 

3 

^ Pm ' ■ 

(31) 


X = 


Y = 
Z = 


H ’ 
H 


Q=Ya = q(^Z), (35) 

N = In a , (36) 

(37) 


With the aid of these definitions, we can write the 
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equations of motion in the autonomous form with ar¬ 
bitrary potentials: 


with 


dX 

In 

In 

dZ 

dN 


= -3X-Y^Q-X^X 

= 


(38) 


^ = -3X2-^(l+a;™)(l-X2-r2) . (39) 


We note that Q is the function of Z. Therefore, Eqs. (1551) 
is indeed an autonomous system of equations. In Table 
II, we present the properties of the five fixed points for the 
quintessence. The points (a, c, e) correspond to the false 
vacua (cti, a^) and real vacuum (< 73 ). The three points 
are stable spirals. On these epoches, the quintessence 
behaves as a damping oscillator with the equation of state 
of firstly behaving as the dust matter and then oscillating 
approaching —1. The points (b, d,) correspond to the two 
local maximum (^ 1 , ^ 2 ) and they are saddle points. On 
these points, the speed of the scalar field exactly vanishes 
and the quintessence acquires the equation of state of 
cosmological constant. 

In Fig. (7-9), we plot the evolution of the rescaled 
speed, X of the quintessence with The figures 
show that, with the increasing of initial speed, the 
quintessence is frozen in the first false vacuum, the second 
false vacuum and the real vacuum, respectively. Com¬ 
pared to evolution of DBI field, we see from Figs.(2-4) 
that the DBI field has more times of oscillations than 
quintessence. How to understand this point? The equa¬ 
tions of motion tell us the friction term due to Hubble 
expansion is 3i70 and 3i70 for quintessence and 

DBI field, respectively. Then with the increasing of speed 
(/), the friction term of DBI field is greatly decreased. So 
the DBI field acquires more times of oscillations. 


IV. CONCLUSION AND DISCUSSION 



FIG. 6 : The quintessence potential V(<?!>) with respect to 
(j>. There are two local maximum (^1 = 0.813, ^2 = 
1.862), one real vacuum ( 0-3 = 2.342) and two false vacuum 
(cti = 0.545, CT 2 = 1.355) for the potential. Physically, the 
quintessence would firstly damped oscillates between these 
vacua and finally dwell on one of the vacuum due to the Hub¬ 
ble friction. 



In general, the equations of motion for generalized DBI 
field and quintessence are rather complicated. So one re¬ 
sort to the method of phase analysis by writing the equa¬ 
tions of motion in the autonomous form. Many special 
form of potentials have been studied for DBI field [ 2 ^ 
and quintessence [T2l - [l3 . [l9l - [^ . However, the general 
method for dealing with arbitrary potentials have not 
yet been proposed. Thus the outcome of this article is 
that we have found the method. 

Different from the potentials studied in Refs. [H- 
[13, [igl - l^ . we investigate the cosmic evolution of the 
generalized DBI field and quintessence with the potential 
of multiple vacua. We find that, with the increasing of 


FIG. 7: Evolution of the rescaled speed {X) of quintessence 
with (j). The point ((/> = cri, V = 0) is a stable spiral and thus 
an attractor. The quintessence is finally frozen at (p = ai. 


initial speed, both generalized DBI field and quintessence 
are successively frozen, in the first false vacuum, the 
second false vacuum and the real vacuum, respectively. 
Compared to the evolution of quintessence, the general¬ 
ized DBI field has more times of oscillations. The reason 
for this point is that, with the increasing of speed 0 , the 
friction term of generalized DBI field is greatly decreased. 
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[Name 

|X 

|y 

^1 <l> 1 

Stability ] 


|W0| 

1 (a) 

1 0 

j 1 

0.545 jcTi = 0.545 [Stable spiral (attractor) j 

1 

i-li 

1 (b) 

1 0 

j 1 

0.813] Cl = 0.813] 

Saddle point ] 

1 

i-li 

1 (c) 

1 b 

j 1 

1.355] (72 = 1.355 [Stable spiral (attractor)] 

1 

i-li 

1 (d) 

1 0 

j 1 

1.862 ]C2 = 1.862] 

Saddle point ] 

1 

i-li 

1 (e) 

0 

1 

2.342 ](T3 = 2.342 [Stable spiral (attractor) ] 

1 

-1 


TABLE II: Properties of the critical points for the quintessence potential given by Eq. (26). 



FIG. 8: Evolution of the rescaled speed (X) of quintessence 
with (j). The point (<() = CT2, X = 0) is a stable spiral and thus 
an attractor. The quintessence is finally frozen at (() = CT2. 


Thus the generalized DBI field acquires more times of os¬ 
cillations. 

The conclusion in this paper may be trivial, but the 
proof and the method are not. As an example, one 
could study the Coleman-De Luccia tunnelling using this 
method. After making the Wick rotation, t = —it 
in Eqs. (21) and Eqs. (38), we are able to study the 
Coleman-De Luccia tunnelling numerically but exactly. 
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